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Abstract. We have studied the magneto-elastic coupling in the double layered Mn perovskite 
Lai.2Sri.gMn207 with an ultrasonic velocity technique. The temperature profile of both the in-plane and 
out-of-plane longitudinal velocities showed a large stiffening anomaly below the insulating paramagnetic 
to metallic ferromagnetic transition. Magnetic fluctuations effects consistent with the layered structure are 
evidenced as a frequency dependent velocity softening above the transition. The magneto-elastic coupling 
has been studied in magnetic field values up to 8 Tesla: the observations are consistent with a substantial 
magnetic anisotropy and a ferromagnetic order parameter with moments lying in the layers. 

PACS. 75.30.Vn Colossal magnetoresistance - 62.65.+k Acoustical properties of solids 

1 Introduction 

Doped perovskite manganites Lai-^A^MnOa (A = Sr, 
Ca) are continuing to attract much attention mainly be- 
cause they are showing a wide variety of magnetic-field 
induced phenomena like the colossal magnetoresistance 
(CMR). In these materials, although the double- exchange 
. mechanism Q| serves as a starting point to explain the 
semiconductor (paramagnetic) to metallic (ferromagnetic) 
phase transition and the associated CMR, the Jahn- Teller 

■ effect and strong electron-phonon interaction must be in- 
| eluded into the model Q . Recent research in these mixed 

valence manganites has expanded to include the search 
for materials with other structure types exhibiting CMR. 
For example, one can adjust the dimensionality of the 

■ Lai-^Sr^MnOa perovskites by inserting a blocking layer 
\ (Lai_ a; Sr a ;)202 between Mn02 sheets to yield the Rud- 

dlesden Popper series of compounds with reduced dimen- 
' sionn, (Lai_ 2; Sr :E )„ + iMn n 03 rl+ i Hf|||,]||. The compounds 
Lai^^Sr^MnOa belong to this family with n = oo. The bi- 
layer perovskite (Lai_ x Sr a; ) n+ iMn207 in which two MnC>6 
layers are stacked with (Lai-zSr^O^ layers along the c 
axis of the structure is obtained for n = 2. The reduced 
dimensionality has been shown to have interesting conse- 
quences on the physical properties of these compounds. 
The material is antiferromagnetic (AFM) for x — 0, while 
it is ferromagnetic (FM) and exhibits a metal insulator 
(MI) transition at T c in the region 0.2 < x < 0.4. The 
maximum T c = 125 K is obtained for x ~ 0.4 and a large 
value of CMR is obtained at temperatures near and far 
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away from T c One of the most distinct features for 

the bilayered manganites is its anisotropic characteristics 
in charge transport and in magnetic properties. They offer 
then a rich opportunity to explore the interplay between 
spin, charge and lattice degrees of freedom in reduced di- 
mensions. In these materials the magnetization process 
causes not only a large negative magnetoresistance but 
also a gigantic anisotropic lattice striction that reveals a 
complex unconventional spin-charge- orbital coupling j^, 
|l^,|lj. This magneto-elastic coupling is dependent on the 
magnetic structure and critical phenomena and its role in 
the mechanism of CMR enhancement has to be clarified. 

Neutron diffraction studies |l2| have revealed that the 
low temperature magnetic phase of these compounds con- 
sists of planar ferromagnetic (FM) and A-type antiferro- 
magnetic (AFM) components, indicating a canted AFM 
ordering. For x — 0.4 the canting angle betwen planes is 
6.3° (nearly planar FM). Hirota et al. also found that 
the A-type AFM ordering remained above T c and that it 
showed an anomalous exponential decrease to Tjv ~ 200 
K. These authors suggested that this AFM phase could 
play a significant role in the enhancement of CMR effects 
in this layered Mn system. In another neutron experiment 
[ p^[ , it was shown that the ferromagnetic transition at 
T c which is exhibiting critical scattering and a divergent 
coherence length is accompanied by 2D ferromagnetic cor- 
relations over a wide temperature range. These FM cor- 
relations were found to fluctuate rapidly close to T c . It 
was then concluded that the layered manganites provide 
a unique opportunity to examine in detail the mechanism 
of the crossover between 2D and ID magnetism. 
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In this paper, we investigate with an ultrasonic veloc- 
ity technique the magneto-elastic coupling in the layered 
compound La1.2Sr1.gM11.2O7 which is showing the maxi- 
mum T c . Because of a strong coupling between spins and 
longitudinal acoustic waves, large anomalies are observed 
on the temperature profile of the acoustic velocity. These 
anomalies have been studied for different crystal directions 
in magnetic field values in the range 0-8 Tesla. These data 
will be discussed in relation with the magnetic structure 
and with the existence of 2D fluctuations on a wide tem- 
perature range. A temperature-magnetic field diagram is 
proposed. 

2 Samples and Experiment 

Single crystals of Lai.2Sri. 8 Mn 2 07 were grown from the 
melt by a conventional floating-zone method using a mir- 
ror furnace E3. The compounds crystallize in the body- 
centered tetragonal structure Ii/mmm (-D4/J JTH]. The 
crystals could easily be cleaved to lead to shiny surfaces in 
the ab plane with the c direction perpendicular to it. Neu- 
tron diffraction experiments on cm size crystals confirmed 
the single crystallinity of the sample. For the present ex- 
periment, the quality of the single crystal, sliced from a 
bigger sample, was also checked with optical polarizing 
microscopy. EDX analysis indicated a homogeneous sam- 
ple without any extra phases (limit of detection 1%) and 
further the composition was found to be close to that of 
the starting material. The ultrasonic velocity is measured 
with a pulsed acoustic interferometer |L6| yielding a sen- 
sitivity in relative velocity variation better than 1 ppm. 
Although the crystals have large dimensions, so large par- 
allel faces that arc well suited for pulsed echo ultrasonic 
measurements, they cannot be used directly. Indeed, the 
magneto-elastic coupling is so strong and the attenuation 
(decibel/mm) in the MHz range so high that we had to 
use much smaller crystals and, doing so, we had to insert 
a delay line (CaF2) to separate the ultrasonic echoes. The 
crystal used in our experiment had dimensions 0.95, 1.15 
and 0.75 mm respectively along the a, b and c axes. One 
parallel face of the crystal is first glued on the surface of a 
CaF 2 delay line (buffer length ~ 7 mm) . Then, a LiNb03 
piezoelectric transducer bonded on the other parallel face 
generates longitudinal waves, at 30 MHz and odd over- 
tones, that will propagate through the crystal-delay line 
ensemble and will be detected by a second piezoelectric 
transducer at the other end of the buffer. The longitudi- 
nal waves can be propagated either along the a (or b) axis 
or the perpendicular direction c. Transverse waves have 
not been used for this study since they could not be iso- 
lated easily from other modes in the ultrasonic experiment 
because of the small dimensions of the crystal. 

The temperature profile of the velocity is obtained by 
monitoring the phase (</>) of the transmitted signal as a 
function of the temperature. However, since two terms, 
<f> = kilx + k^li where hi, k%, l\ and I2 are respectively 
the ultrasonic wave vectors and the lengths of the crystal 
(1) and the delay line (2), contribute to the total phase 
of the transmitted signal, it is necessary to substract the 
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Fig. 1. Relative variation of the longitudinal velocity as a 
function of temperature obtained at 165 MHz: crystal/buffer 
(full line) and buffer (dotted line). 

delay line contribution fah) to isolate the crystal one 
(kxli). The former is measured in a pulsed reflection ex- 
periment in the delay line only at the same frequency 
and for identical experimental conditions. Our ultrasonic 
data have not been corrected to take into account ther- 
mal expansion (striction §,0) effects (^) since these 
are generally orders of magnitude too small compared to 
the measured relative velocity variation (— ). The abso- 
lute values of the longitudinal velocities are around v Qj t> 
= 4300 and v c = 3400 m/sec. No anisotropy was found 
between the a and b directions. The temperature profile of 
the velocities is obtained in the range 4-200 K in magnetic 
field values between 0-8 Tesla. The temperature is moni- 
tored either with a Si diode or a carbon glass sensor and 
stabilized with a LakeShore controller. For reasons of ex- 
cellent reproducibility during thermal cycling, silicon seal 
was used for all bonds between the crystal, the transduc- 
ers and the buffer; no velocity data can be obtained above 
200 K where the silicon seal presents a phase transition. 



3 Results and Discussion 

In figure 1 we present an example of the temperature pro- 
file of the variation of the velocity (relative to its value at 
200 K) for 165 MHz longitudinal waves propagating either 
trough the crystal/buffer ensemble (full line) or the buffer 
alone (dotted line). The velocity in the buffer shows the 
usual monotonous behavior which is related to the an- 
harmonic contribution. This means that all the features 
observed on the crystal/buffer profile are due to magneto- 
elastic effects in the Lai.2Sri. 8 Mn 2 07 crystal. When the 
buffer contribution (fc2^) to the overall phase (0) is sub- 
stractcd, wc get the temperature profile of the crystal ve- 
locity that is shown in figure 2 for 165 MHz longitudinal 
waves having their wave vector k parallel to the a (or b) 
axis or perpendicular to it (k || c). When the temperature 
is decreased below 200 K for both crystal directions, the 
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Fig. 2. Relative variation of the longitudinal velocity as a 
function of temperature obtained at 165 MHz: k || a or b (up- 
per panel) and k \\ c (lower panel). The arrows indicate the 
transition temperature T c . 
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Fig. 3. Microwave resisitivity (17 GHz) along c for the layered 
manganite Lai.2Sri.gMn207 as a function of temperature. The 
dashed line indicates the transition temperature T c . 



obtained at 17 GHz on the same crystal. The semicon- 
ductor metal transition is characterized by an abrupt de- 
crease of the resistivity at T c = 125 K where a dashed 
line indicates the maximum variation of the resistivity. 
The temperature profile and the transition temperature 
are similar to the DC ones reported in the literature |^|. 
In the ultrasonic experiment, this value of T c corresponds 
to the temperature of maximum stiffening rate observed 
below the weak softening peak as indicated by an arrow 
in Fig. 2. 



velocity first increases smoothly but, in the vicinity of the 
phase transition temperature T c , a small softening of the 
velocity is observed before an important stiffening progres- 
sively sets in below T c . At low temperatures the smooth 
behavior is recovered although a small slope increase is 
seen around 20 K. Differences between the two directions 
are a larger stiffening anomaly below T c and a more pro- 
nounced softening around T c when waves propagate along 
the a (or £>)axis. As it is observed on the transport proper- 
ties Q , when a 8 Tesla magnetic field is applied along a (or 
b), the stiffening anomaly is shifted to higher temperatures 
as shown in Fig. 2. Such a large acoustic anomaly (~ 2% in 
amplitude) is typical of what is found in strong magnetic 
materials and we thus suggest to relate it to a magneto- 
elastic coupling at the paramagnetic-ferromagnetic transi- 
tion. Indeed, since quasi-particle screening effects are ex- 
pected to be orders of magnitude smaller and to rather 
yield a softening anomaly when entering the metallic state 
below T c , the charge degrees of freedom cannot explain 
the features observed here. 

In order to interpret the velocity data presented in 
Fig. 2, we need to define correctly the transition tempera- 
ture T c . This can be done most easily on the temperature 
profile of the resistivity p(T) where the maximum rate 
(dp/dT) generally defines T c . We present in figure 3 the 
microwave resistivity along c as a function of temperature 



For T<T C the stiffening anomaly is likely due to a cou- 
pling between ferromagnetic spins and longitudinal acous- 
tic phonons. Indeed, the stiffening anomalies in zero field 
shown in Fig. 2 are coherent with this picture: a sponta- 
neous magnetization progressively builds up below T c and 
this gives an increase of Av/v. We have tried to isolate the 
magneto-elastic coupling by substracting the anharmonic 
contribution which is found by extrapolating to low tem- 
peratures the profile observed for T>T C . The result is pre- 
sented in figure 4 for both crystal directions. A stiffening 
anomaly which mimics the magnetization M of the ferro- 
magnetic phase is obtained. At T c , the stiffening appears 
to be sharper along the c axis than along a (or 6) and the 
anomalies have different amplitudes, being larger along 
the ab plane. A change of slope below 20 K is clearly ob- 
served on both curves. Its origin could not be identified as 
it is hardly modified by a magnetic field and, thus, it will 
not be discussed further. The softening observed above T c 
is due to magnetic fluctuations and this will be discussed 
shortly. 

The stiffening anomaly can be explained by a magneto- 
elastic coupling according to the following scenario. We 
consider a coupling between the strain e and the order pa- 
rameter Q (spin-spin correlation function) and the energy 
term is phenomenologically expanded in powers of these 
parameters [^) : a stiffening anomaly below T c is obtained 
if we consider a bi-quadratic coupling energy term Xe 2 Q 2 . 
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Fig. 4. Magneto-elastic contribution to the velocity profile at 
165 MHz in zero field: a or & axis (full line) and c axis (dotted 
line). Inset: c axis data as a funtion of the reduced temperature 
(1 - £-) on a log- log scale; the straight line represents the best 
fit in the critical region with j3 = 0.12. 



In most favorable cases where a single ferromagnetic com- 
ponent is present in the low temperature phase and where 
domain effects are absent, the measurement of the elastic 
constant C, so the velocity v (C = pv 2 , p being the mass 
density), can reveal the temperature dependence of the 
order parameter. The curves shown in Fig. 4 are in agree- 
ment with such a scenario. Since the largest anomaly is ob- 
served for in-plane propagation, a larger coupling constant 
A is expected. However, ferromagnetic domains are likely 
present in Lai^Sri.gM^Or crystals and a AFM compo- 
nent has also been identified in neutron experiments ||12f 
at low temperatures; we thus have to be careful to extract 
the order parameter from these curves, especially in the 
critical temperature region. Indeed, the data presented in 
figure 2 reveal that a 8 Tesla magnetic field oriented along 
a (or b) modifies differently the low temperature ampli- 
tude of the anomaly according as the wave is propagating 
parallel or perpendicular to the plane. For the parallel 
case, a further stiffening likely due to domain wall effects 
is observed in the ferromagnetic phase when no effects are 
detected for the perpendicular one. 

Longitudinal waves can interact most easily with do- 
main walls that are perpendicular to the elastic polariza- 
tion. For a tetragonal structure where the Mn magnetic 
moments are ordered ferromagnetically along the a or b 
axes, four domains are expected in zero field and the walls 
should run preferentially along the c axis. The application 
of a magnetic field will affect the domain structure and 
the velocity should be changed accordingly. These mag- 
netic field effects investigated at 10 K, far below the crit- 
ical region, are presented in figure 5 for different configu- 
rations of wave propagation and field orientation. In the 
upper panel we are considering waves propagating along 
a (or b). When the field is along the same direction (easy 
plane), the velocity first softens weakly, reaches a mini- 
mum around 0.12 Tesla, then stiffens and saturates above 
0.4 Tesla. The minimum could be the result of an increase 
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Fig. 5. Variation of the velocity relative to its value at zero 
field at 165 MHz as a function of magnetic field at 10 K. Upper 
panel, wave vector k || a or b; lower panel, k || c. 



softening due to domain wall motion before stiffening is 
obtained when the 90 degrees domains rotate to align 
along the field. When the field is rather along the hard 
axis (c), no minimum is observed and saturation is ob- 
tained at a much higher field value around 1.3 Tesla. This 
is consistent with a large magnetocrystalline anisotropy 
with the hard axis along the c direction; as the domain 
structure is not modified by a perpendicular field, no soft- 
ening is observed for this configuration. In the lower panel 
of Fig. 5, we notice that longitudinal waves propagating 
along c are not affected if the field is along a (or b) since 
interaction with domain walls is not efficient. For field 
along c however, softening is obtained when the moments 
are forced to leave the easy plane with a maximum ob- 
tained around 1.4 Tesla and saturation is observed above 
2 Tesla when all the moments are aligned in the field. The 
negative value obtained at saturation is likely related to 
the demagnetization field. This dependence of the velocity 
in the ferromagnetic phase is fully coherent with the mag- 
netic structure and the field behavior observed by others 
|l8| although the absolute values are somewhat different. 
In these measurements however, the possible role played 
by a small AF component is difficult to evaluate. 

The effect of domains on the velocity data complicates 
the determination of the temperature profile of the order 
parameter, especially the critical behavior near the phase 
transition temperature. If a bi-quadratic magneto-elastic 
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Fig. 6. Relative variation of the velocity as a function of tem- 
perature near T c at different frequencies in the MHz range and 
in zero field. Upper panel, wave vector k || c; lower panel, k \\ a 
or b. 



energy term is used, one can show that the relative velocity 
chan ge i s proportional to the square of the order param- 
eter [p7| . Now, if we consider the c axis velocity data of 
Fig. 4 as being free of magnetic domains effects, the critical 
region near the phase transition can be fitted reasonably 
well to the following relation, (^)s ~ Q(T) ~ (T c -T)P, 
with parameters T c = 125.5 ± 0.2 K and (3 — 0.12 ± 0.01 
(see inset of Fig. 4). The small value of the exponent is 
consistent with the presence of strongly two dimensional 
fluctuations below T c as predicted by the 2D Ising model 
(P = 0.125). If this weak value of (3 agrees with the neu- 
tron scattering data of Osborn et al. |13] and Rosenkranz 
et al. pp[| , it contradicts the 3D Heisenberg value (0.35) of 
Chatterji et al. However, our results indicate clearly 
that 2D magnetic fluctuations are dominating the critical 
temperature region. At much lower temperatures, the c 
axis data of Fig. 2 (lower panel) reveal that, although a 
magnetic field shifts the phase transition to higher tem- 
peratures, it has no effect on the low temperature value 
of the order parameter Q(0). Such a field independent or- 
der parameter Q(0) has also been observed previously in 
low-dimensional AFM systems p2],p3[. 

The presence of 2D magnetic fluctuations can also be 
inferred directly from the velocity data presented in figure 
4. Indeed, the velocity softening observed just above T c 
for in-plane propagation can be attributed to fluctuations 
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Fig. 7. Relative variation of the velocity at 165 MHz as a 
function of temperature near T c for different magnetic field 
values (H |j a or b): wave vector k oriented along a or 6 axes. 



and the fact that this softening is hardly seen for c axis 
propagation implies that the fluctuations possess a strong 
2D character. In fact, the fluctuation regime for T >T C is 
highly dependent on the frequency in the MHz range. This 
can be noticed in figure 6 where we present the velocity 
softening obtained at different ultrasonic frequencies near 
the transition. For in-plane propagation (lower panel), the 
amplitude of the softening increases dramatically when 
the frequency is decreased from 165 to 32 MHz (it reaches 
more than 2% at the lowest frequency, the same value 
as the stiffening due to the order parameter below T c ) 
and its maximum comes very close to T c (125 K). For 
out-of-plane propagation (upper panel), a similar effect 
is obtained but the softening amplitude is much smaller 
by at least one order of magnitude. This clearly indicates 
that, for T >T C , 2D fluctuations are strongly coupled to 
acoustic phonons and that their relaxation time is larger 
than 10~ 7 sec. The range of temperatures over which these 
fluctuations are extending is difficult to determine since 
we do not know exactly the temperature profile of the 
velocity without the effects of spins. According to the data 
of Fig. 6 these fluctuations can be detected at least up to 
150 K. This frequency dependence is consistent with the 
ferromagnetic character of the 2D fluctuations since these 
are expected to be enhanced as the wave vector q — > 0. 

It is also well known that the application of a magnetic 
field can accentuate the fluctuations in low-dimensional 
magnetic systems. This can be seen in figure 7 where we 
indicate how the critical region for in-plane propagation is 
modified by a magnetic field oriented in the plane. In fact, 
the field is not only shifting the transition temperature to 
higher values, but it also increases its width. Besides, the 
amplitude of the softening increases first with field up to 
1 Tesla and then decreases smoothly as the critical re- 
gion becomes wider and wider. The same trend is also 
observed for out-of-plane propagation with much smaller 
effects. This behavior is again consistent with the pres- 
ence of 2D ferromagnetic fluctuations strongly coupled to 
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Fig. 8. Magnetic phase diagram for both magnetic field orien- 
tations: T C (H) is determined from the maximum slope of the 
velocity variation in the critical region for both wave propaga- 
tion configurations. 



acoustic phonons over a wide temperature range above 
T c . Whether the character of the fluctuations is 2D Ising 
is not fully established since an effective finite-size 2D XY 
model has been suggested to explain quantitatively the 
critical properties above and below T c |pof i although the 
exponent predicted by this model seems too large (0.23) 



Finally, the shift of the transition temperature T c as a 
function of the magnetic field has also been investigated 
for both wave propagation configurations. The deduced 
magnetic phase diagram is shown in figure 8. The critical 
temperature T C (H) has been determined with the crite- 
rion used in zero magnetic field. As expected, there is es- 
sentially no difference between the T c 's determined from 
in-plane and out-of-plane propagation for both orienta- 
tions of the magnetic field, but only for H < 3-4 Tesla. 
For higher field values, the T c 's differ surprisingly. We will 
thus discuss these field ranges separately. When < H < 
3 Tesla, the magnetic phase diagram is not dependent on 
the choice of the wave propagation but it clearly depends 
on the field orientation. When the field is oriented along 
one of the easy axes a or b (lower panel), T C (H) increases 
quasi- linearly with field up to 1.5 - 2 Tesla with a rate 
of ~ 8 K/Tesla; for higher fields, the rate gets smaller 
around 6.3 K/Tesla. For the field oriented along the hard 
axis c (upper panel), the relation between T c and H is 



rather quadratic at low field values as indicated by the 
positive curvature noticed in Fig. 8. At higher fields (H > 
1.5 Tesla), we get a linear increase with a rate similar to 
the easy plane orientation one. This slower increase of T c 
along the c axis is coherent with an anisotropy field around 
1.5 Tesla as determined previously from the field depen- 
dence of the low temperature stiffening anomaly (Fig. 5). 
For 3 < H < 8 Tesla, if T c is determined from the c axis 
propagation, the linear dependence on field observed be- 
low 3 Tesla is maintained up to the maximum value of 8 
Tesla. T C (H) determined from in-plane wave propagation 
is however smaller. Since the critical region is much less 
affected by the fluctuations and domain walls for wave 
propagation along the hard axis, the true T C (H) must be 
determined from this configuration; in the case of wave 
propagation along one of the easy axis, fluctuations widen 
the transition in such a way that T c appears to be under- 
estimated. The open circles of figure 8 are thus fully rep- 
resentative of the true phase diagram for this compound. 



4 Conclusion 

We have identified in the layered Lai^Sri.gM^Cv CMR 
compound an important coupling between longitudinal 
acoustic phonons and ferromagnetic moments. In the fer- 
romagnetic phase at low temperatures, this coupling pro- 
duces a stiffening anomaly on the ultrasonic velocity from 
which the temperature profile of the order parameter, the 
spin-spin correlation function, could be determined if do- 
main walls effects are taken into account. The field de- 
pendence of this stiffening anomaly at low temperatures 
confirms the known magnetic structure of this compound: 
an easy plane of magnetization along ab and a hard axis 
along c with an anisotropy field just above 1 Tesla. In the 
critical region, the exponent of the temperature power law 
indicates that a 2D Ising model is appropriate to describe 
the magnetization process. A small anomaly was also ob- 
served for T< 20 K but its origin could not be identified. 
We have also observed above T c an important softening of 
the velocity which confirms that 2D ferromagnetic fluctu- 
ations are present in the paramagnetic state over a wide 
temperature range. However, with our ultrasonic measure- 
ments, we could not identify any effects due to the pres- 
ence of an antiferromagnetic order or fluctuations. 
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